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ABSTRACT 

Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the 
high-density regions within molecular clouds. 

Aims. We measure the mid-infrared extinction law produced by dense material in molecular cloud cores. Since the extinction at these 
wavelengths is caused by dust, the extinction law in cores should depart from that found in low-density environments if the dust grains 
have different properties. 

Methods. We use the unbiased LINES method to measure the slope of the reddening vectors in color-color diagrams. We derive the 
mid-infrared extinction law toward the dense cores B59 and FeSt 1-457 in the Pipe Nebula over a range of visual extinction between 
10 and 50 magnitudes, using a combination of Spitzer/IRAC, and ESO NTT/VLT data. 

Results. The mid-infrared extinction law in both cores departs significantly from a power-law between 3.6 and 8 /urn, suggesting that 
these cores contain dust with a considerable fraction of large dust grains. We find no evidence for a dependence of the extinction law 
with column density up to 50 magnitudes of visual extinction in these cores, and no evidence for a variation between our result and 
those for other clouds at lower column densities reported elsewhere in the literature. This suggests that either large grains are present 
even in low column density regions, or that the existing dust models need to be revised at mid-infrared wavelengths. We find a small 
but significant difference in the extinction law of the two cores, that we tentatively associate with the onset of star formation in B59. 

Key words. Radiation mechanisms: general; Methods: observational; Techniques: photometric; Stars: formation; ISM: clouds; 
Infrared: ISM 



1. Introduction 

Dust plays a crucial role in many of the chemical and physical 
processes that take place in the interstellar medium. It plays an 
especially important role in molecular clouds and in the forma- 
tion of stars. At low densities, dust grains are the catalysts for 
the formation of molecules that then cool the clouds radiatively 
(e.g., |Whitworth et al.|1998||Draine|20Tl) . At high densities the 
molecules begin to deplete from the gas phase and dust becomes 
the dominant coolant, lowering the gas temperature enough that 
fragmentation of cores into stars may occur ( W hitworth et al.| 
1 1998 ; Larson 2005). Yet, the properties of the dust grains, in- 
cluding their composition and size distribution, and how they 
evolve in different environments are not accurately known. 

As the density increases, dust grains are believed to grow by 
coalescence through grain-grain inelastic collisions and by the 
accretion of ice mantles, as molecules freeze out on to their sur- 
faces (e.g., |Whittet et al.|1988||Ossenkopf|1993| l. The ice man- 
tles are typically very thin ( |Draine|1985| l and so do not increase 
the grain size significantly by themselves, but they promote fur- 
ther growth by making the grains "stickier" in collisions (e.g., 
Orm el et al.|200 9). The size distribution of dust grains is there- 
fore expected to vary from the low-density interstellar medium 
to the cold, dense environment of cores as the small grains and 



* Based on observations made with ESO Telescopes at the La Silla 
Paranal Observatory under programme IDs 62.1-0866, 65.1-0083, 67. C- 
0514, and 069.C-0426. 



PAHs (polycyclic aromatic hydrocarbons) disappear and popu- 
late the larger grain-size side of the distribution. 

The impact of grain growth on star formation is not clear, but 
there are indications that it may play a part, 
gests that the grain size distribution is not rel 
of temperature at the highest densities (> 10~ 



Larson (2005 ) sug- 



evant to the control 

-18 



g cm ), but that 
it can be important at intermediate densities where it can still af- 
fect the collision rate between the gas and dust. This defines the 
density at which gas and dust become thermally coupled, which 
according to the author ultimately determines the characteristic 
stellar mass. |Omukai (2000), on the other hand, suggests that 
different grain size distributions can produce different minimum 
temperatures, again influencing the fragmentation properties of 
the cloud, although, according to Larson ( |2005| l, their models 
produce unrealistically low temperatures considering the radia- 
tive heating by ambient far-infrared radiation. 

The extinction law is often used as a proxy for the size dis- 
tribution and composition of dust grains, traditionally in the op 
tical and ultra-violet (e.g., |Savage & Mathis|1979|[Draine|201 1 



and more recently also in the near- and mid-infrared. The op- 
tical and ultraviolet have been successfully used for decades to 
study the low-density interstellar medium, revealing an extinc- 
tion law that can be described as a function of a single param- 
eter, the total- to-selective extinction ratio, Ry = Av/E(B - V) 
dCardelli et aT1[T989) |Qlofsson & Olofsson||20T0l i. In the near- 
infrared (NIR) the extinction law is considered invariant and 
well described by a power law {A^ AT") of index a in the 
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range 1.6- 1.85 (]Koomneef|1982||Wh7ttet et al.|1988[|Car delli 



et al. 



1989 



Martin & Whittet||1990| |Clayton & Mafhisl|1988 
Becklin et al.||19781|Landini et al.||19841 [Rosenthal et al.pOOO 
Nishiyam a et al.||2006[ |Olofsson & Olofsson| |2010| l, and in the 
mid-infrared (MIR, 3-8 jum in this paper), the wavelength range 
of interest for the present study, the extinction law is not yet 
well-characterized or understood. Steinack er et al. | ((20 1 ) have 
proposed another tool to probe the grain size distribution of dust 
using scattering of light rather than absorption. They used the 
"coreshine" method to infer the presence of large (up to 1 fim) 
grains in dense molecular cores. 

In this paper we compare the infrared extinction law in two 
dense cores of the Pipe Nebula: FeSt 1-457, which is quiescent, 



and B59, which contains a small young cluster (Forbrich et al. 
2009; Covey et al.|20 10l. The Pipe Nebula is one of the closest 



molecular clouds to the Sun (d — 130 pc, Lomba rdi et al.|20 06a ) 
and has been studied extensively over the few past years. The 
first systematic study of the cloud belongs to Onishi et al. ( 1999| l, 



who present CO observations of a large part of the clo ud and esti- 



mate the total 12 CO mass to be of the order of 10 4 ML. Lombardi 



resolution extinction map based on 2MASS data covering 8°x6° 
later improved by Roman-Zuniga et al. (2010) using deep in 
frared data 



Alves et al. (2007 



identify and characterize 159 
cores within the Pipe based on the extinction map, ranging from 
3.0 to 24.3 visual magnitudes of extinction, and with correspond- 
ing masses between 0.5 to 28 Mq, and use this information to 
derive the Pipe's dense core mass function. Mue nch et al.| (f2007 ) 
establish most of these cores as actual dense structures, rather 
than the superposition of less dense material, using observations 
of the 7=1-0 transition of C ls O, further analyzing the cores's 
kinematical structure. Rathborne et al.| ([2008 ) further validate 
the association of the cores with dense gas using molecular line 
observations of NH 3 (1,1), NH 3 (2,2), CCS (2i-l ), and HC 5 N 
(9,8). Using the C ls O radio data and the data from the extinction 



map, Rathborne et al. (2009J produce a refined sample of cores 
and confirm the previously found core mass function. Frau et al. 
(2010| and Frau et al. ( 2012| l study the size, density and chem- 
istry of a few cores in the Pipe using radio data, proposing an 
evolutionary scenario for the chemistry in the cloud. The mag- 
netic fields of the Pipe were studied by Alves et al. ( 2008 ) using 
optical polarimetry. They divide the cloud in three zones accord- 
ing to the properties of the magnetic field: B59, the densest core 
in the cloud, the long "stem", and finally the "bowl", with the 
strongest magnetic field (15% polarization). They estimate the 
mass-to-flux ratio to be approximately super-critical toward B59 
and sub-critical inside the bowl. Finally, Peretto et al. ( 2012[ i and 
Gritschne der & Lift] ( |2012[ ) address the formation of the cloud 
based on Herschel data and analytical calculations. The extinc- 



tion law in the Pipe was first studied in the NIR by Lombardi 



et al. (2006b i, who find an average normal extinction law for the 
whole cloud using 2MASS data. Roman-Zuniga et al. (2007) 



study the extinction law of the B59 core in particular and find it 
to be flat in the 3.6-8.0 yum regime, while confirming a normal 
law in the NIR. 

In this paper we will readdress the extinction law of B59 us- 
ing the LINES method presented in |Ascenso et al. (2012 here- 
after referred to as Paper I), and determine the extinction law 
of another dense core in the Pipe, FeSt 1-457. This core is lo- 
cated in the opposite extreme of the nebula from B59, and has 
no signs of active star formation. Agu ti et al.| ( |2007| ) suggest the 
core is gravitationally bound in a quasi-stable state near hydro- 
static equilibrium, while the outer layers of the core are proposed 
to be oscillating. 



Studying the extinction law in the context of the dust prop- 
erties in two cores of the same cloud has the clear advantage 
that the chemical composition and in general the environment 
and conditions external to the cloud must be the same. The dif- 
ferences between the two cores can therefore be isolated and 
studied in detail. 

In Section [2] we present the data used for this study, and in 
Section[3]we derive the extinction laws. In Section|4]we discuss 
the results in the context of previous studies and models, and we 
present our conclusions in Section [5] 

2. Observations 

This section describes the data used to derive the extinction 
law in B59 (J2000, 17hllm21s, -27°26'10") and FeSt 1-457 
(J2000, 17h35m47s -25°33'09"), in the Pipe Nebula molecular 
cloud. The data used here will be made available in the context 
of a forthcoming paper. 



et al.|([2006a[) refine the distance to the cloud and present a high- 2. 1 . B59 



We used the dataset for B59 of Roman-Zuniga et al. ( 2007 ). The 
near-infrared data were obtained with ESO NTT/SOFI and ESO 
VLT/ISAAC in the H (1.6 jum) and K s (2.2 /zm) bands. The / 
band (1.2 /im) was only partially observed, and the catalog was 
therefore complemented with 2MASSQdata for better coverage. 
This added 162 sources with respect to the original catalog of 
Roman-Zuniga et al. (2007). Unfortunately, because 2MASS is 



too shallow, there are still no sources with J band information in 
the center of the core. 

The mid-infrared (Spitzer/IRAC) dataset was the exact same 
as inll 



i Roman-Zuniga et al. (2007): the publicly available catalog, 
pan of the Cores to Disks Legacy Project ( Evans et al.|2003| ). 



2.2. FeSt 1-457 
2.2.1. NIR data 

The main observations for FeSt 1-457 were made with the 
ISAAC and the SOFI near-infrared imagers, mounted, respec- 
tively, at the 8.2 m UT3 unit of the ESO/VLT, and at the 3.5 
m ESO/NTT. The FeSt 1-457 observing runs were completed 
in March 1999 (SOFI, program ID 62.1-0866) and July 2000 
(ISAAC, program ID 65.1-0083). All the data are currently avail- 
able in the ESO archives. SOFI has a field of view (FOV) of 
5' x 5' with plate scale of 0.288 arcsec pixeL 1 , almost seven 
times the angular resolution of 2MASS, and higher sensitiv- 
ity: achieved limits are estimated to be about 4-5 mag deeper 
than 2MASS in each band. Such characteristics can resolve the 
densely crowded Galactic Bulge field in the background of the 
Pipe Nebula and penetrate in regions with up to about 50 mag of 
visual extinction. The field of view of SOFI is wide enough to 
completely contain some of the larger dense regions in the Pipe 
Nebula (R ~ 0. 10. 15 pc), with the exception of the central core in 
Barnard 59. The ISAAC observations were intended to fully re- 
solve the centers of the most dense and obscured regions. These 
observations have a resolution of 0.144 arcsec pixeL 1 with a 
field of view a quarter of the size of that of SOFI. 



The Two Micron All Sky Survey is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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2.2.2. MIR data 

The Sjpifeer/IRAcQdataset for FeSt 1-457 is from the Cycle I 
GTO observations and from the Cycle II GO Pipe Cloud Survey. 
The data were available in the archive as processed by the Spitzer 
Science Center team using pipeline version S18.7 that already 
includes artifact correction. Each field was observed with long 
(10.4s) and short (0.4s) exposure times. The source detection 
and photometry were done using the MOPE^] pipeline soft- 
ware APEX multiframe on the artifact-corrected images. The 
main difficulty in this step was the extreme crowding caused by 
the position of the Pipe Nebula against the Galactic Bulge. The 
crowding was severe in all fields except toward the very center 
of the core where the extinction eliminates most of the (back- 
ground) sources. Apart from the technical problems that arise 
when the software attempts to fit several PRFs (Point Response 
Functions) to many "clustered" sources at once, the crowding 
ultimately produced uncertainties in the photometry of sources 
with bright neighbors. Unlike the typical photometric error, esti- 
mated by the pipeline during the detection and fitting procedure, 
these errors are difficult to quantify on an individual basis. Since 
the crowding is a function of wavelength, Channel 1 (3.6 pm) 
being the most crowded due to its higher sensitivity, these un- 
certainties will be larger for the shorter wavelength bands. The 
source extraction and photometry were done independently for 
the long and short exposures. 

Only the objects with the best quality photometry were se- 
lected. The quality criteria were a relative flux error below 10%, 
and ax 2 of the PRF fit per degree of freedom below 25, which 
ensured a clean sample. After the error cuts, the long and short 
exposure catalogs for the same channel were merged such that, 
for the stars common to both, those with the smaller SNR were 
discarded. The final catalogs for the different channels contain- 
ing only the sources obeying these criteria were then merged 
with the NIR data. 



2.3. Control field 

No dedicated control field was taken for this study. Instead we 
used an extinction-free nearby field close to the core B72, the 
same that was used by Roman-Zuniga et al. ( 2QQ7\ in their 
work. For this field the near-IR observations were from the same 



ESO/SOFI survey as described in section 2.2.1 and the mid- 
infrared data were from the Spitzer/IRAC catalog of the Cores 
to Disk (c2d) Legacy project (Eva ns et al.|20 03). The measured 
photometric dispersion for background sources is comparable to 
the uncertainties in the (A - K s ) colors. This indicates that the in- 
trinsic color dispersion in the Galactic Bulge background field of 
the Pipe Nebula is small enough that our photometric accuracy 
cannot resolve it. Because these data were not meant as control 
fields for the present study, they are shallower than the science 
observations. 
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Fig. 1. (K s - A) vs. (H - K s ) color-color diagrams for B59, for 
A = 3.6 pm (top left), A = 4.5 pm (top right), A = 5.8 pm 
(bottom left), and A — 8.0 pm (bottom right). The solid lines 
illustrate the reddening vectors for each band: the slopes are the 
LINES estimates, and the intercepta are such that the lines pass 
through the median of the distributions. 
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Fig. 2. Same as Figure[T|but for FeSt. 



3. The extinction law in B59 and FeSt 

3.1. LINES determination ofA A /A Ks 

The data described above was used to construct (K s - A) vs. 
(H-K s ) color-color diagrams, where A is 3.6, 4.5, 5.0 and 8.0 pm 

2 This work is based in part on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. 

3 http : //irsa . ipac . caltech . edu/data/SPITZER/docs/ 
dataanalysistools/tools/mopex/ 



for the IRAC bands. Since the observed stars are mostly giants 
behind the FeSt and B59 cores, their position in the color-color 
diagrams is dictated by the loci of their intrinsic colors, by their 
photometric error, that scatters the points randomly in the dia- 
gram, and by the extinction caused by the cores, that moves the 
points along the reddening vector. The latter is the dominant ef- 
fect, as can be seen in Figures [T] and [2] and it is the slope yS of 
these vectors in the different color-color diagrams that we wish 
to measure in order to derive the extinction laws. LINES (Linear 
regression with Errors and Scatter) has been validated in Paper I 
as the best and most unbiased method to do this. 
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The LINES method determines the slope of a linear distribu- 
tion of points with errors in both directions, possibly correlated, 
and with an intrinsic distribution that may not itself be linear. In 
the case of the reddening vector in a color-color diagram, the in- 
trinsic scatter is the locus of the stellar population that is being 
reddened by the cloud, and we wish to disentangle this scatter 
from the pure contribution of the extinction. The intrinsic dis- 
tribution of the points is sampled by a control field expected to 



sample the same stellar population. As described in section 2.3 



our control field was observed with a different setup from the 
science field, producing a shallower dataset. However, since the 
background population of the Pipe Nebula is dominated by gi- 
ant stars which have a very narrow range in infrared colors, the 
control field, albeit shallow, represents the intrinsic scatter in the 
science field appropriately. 

The slope fix (- E(K S -A)/E{H-K S )) of the reddening vector 
is calculated as: 



Cov(x,y) - Cov(e x , e y ) - Cov(x cf ,y cf ) + Cov(e cfx , e cfy ) 
Var(x) - Var(e-0 - Var(jc c /) + Var(e cfx ) 



vaiyx) — varie-'j - varix 1 -' j + \w\e^) 

(1) 

where x and y are the (H - K s ) and the (K s - A) colors, respec- 
tively, and Cov(x, y), Cov(e x , e y ), Var(jt), and Var(e JC ) are the co- 
variance of the (x,y) distribution, the covariance of the (x,y) er- 
rors, the variance of the x distribution, and the variance of the 
errors in x, respectively; the symbols marked with cf are the 
same but refer to the control field. The uncertainty in fix is es- 
timated using the bootstrap method corrected by a factor of 1 .25 
as described in detail in Paper I. 

Figures [T] and [2] show the best LINES fits to the color-color 
diagrams of B59 and FeSt, respectively. The fit is naturally 
weighed more heavily toward the low-extinction end, since it 
contains most of the datapoints. The apparent deviation between 
the fit line and the data points at high extinction is not considered 
to be statistically relevant (see Sect. |3.3| >. 

The extinction law Ax/Ak s is derived from fix as: 



Ax_ 

As, 



l 



1A 



(2) 



The uncertainty in the ratios A,\/Ak s is calculated by error prop- 
agation from the uncertainty in fix- We assume a ratio A/y/A^ of 
1 .55 following Indebetouw et al. ( 2005 [ l since we cannot deter- 
mine it accurately from our data. The choice of this anchor point 
is crucial for the comparison between different observations and 
between observations and models, since different values will 
produce different extinction laws. We chose this value over oth- 
ers in the literature (for example, most recently, Nishiyama et al. 
d2009b found A H /A Ks = 1.62 toward the Galactic Center) be- 



cause it is the most widely used in previous studies and therefore 
allows for the most direct comparison. Also, using this value, fi j 
translates into a value of A j/A Ks in very close agreement with the 
models considered here (see Fig. BJ. All these models agree on 
the value of Ay/A^ even though they refer to different grain size 
distributions, so this agreement lends support to our assumption. 
Regardless, our results can be easily adapted to another defini- 
tion should it prove to be more appropriate (see also Sect. |4.3)>. 



3.2. The extinction law in B59 and FeSt 

The derived slopes fi A and the corresponding extinction law 
Ax/Ak s from 1.2 to 8.0 fim are given in Table [T] for B59 and 
FeSt, and graphically in Figure [5] (red and bluelines for B59 
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Fig. 3. The extinction law for FeSt (blue line) and B59 (red line). 
Also shown are the Draine (2003a b) models for Ry = 3.1 (solid 
line) and Ry = 5.5 (dotted line), and the Weingartner & Draine 
( [20011 > models for R v = 5.5, case B (dashed line). 



and FeSt, respectively). Also shown in the figure are the mod- 
els from Weingartner & Draine ( 2001 1 hereafter referred to as 



WD01): the solid and dotted lines represent the "case A" model 
of WD01, re-normalized by |Draine | (2003a bf] for R v of 3.1 



and 5.5, respectively, and the dashed line represents the "case 
B" model of WD01 for R v of 5.5. 

In the near infrared only the J band relative extinction can 
be assessed, since A K /A K = 1 and we fixed A H /A K = 1.55 



(see Sect. 3.1 1. Our derived Aj/A K . for both B59 and FeSt are 



in perfect agreement with the dust grain models, and are also 
indistinguishable from other measurements in the literature for 
this wavelength (see Sect. 4.1 1. This agreement supports our as- 



sumption of Ah/Ak s = 1.55. 

In the IRAC bands, from 3.6 to 8.0 yum, both cores show a 
very flat extinction law, in particular flatter than the Ry = 5.5 
"case B" model ofWDOl. 



3.3. The extinction law as a function of column density 

We analyzed the behavior of the extinction law as a function 
of extinction to see whether we found a change that could be 
attributed to grain growth at any specific column density thresh- 
old. Figures [4] and [5] show that the datapoints at high extinctions 
depart slightly from the LINES fit, suggesting a change of be- 
havior at some column density. In Paper I we showed that the 
LINES method is robust enough to detect such a change. Using 
the (H - K s ) color, we divided the sample in low-extinction 
((H - K) less than a value (H - K)u m i t ), and high-extinction 
((H - K) > (H - K)n m j,) subsamples, and determined the best 
fits to the reddening vector in the two groups using LINES, ob- 
taining two slopes Pi ow and fihigh f° r eacn (H - K)u m i,. This was 
done for increasing values of (H — K)u m i, in steps of 0.1 mag. 
Figures [4] and [5] show the results for B59 and FeSt, respectively. 



4 http : // www . astro . pr inceton . edu/ ~ draine/ dust/ 
dustmix.html 
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Table 1. Slopes E{K S - A)/E(H - K s ) and the extinction law A^/A^ as a function of wavelength for B59 and FeSt. 



B59 FeSt 1-457 



Band 


A (/um) 


fa = E(K S -A)/E(H-K S ) 


A A /A Ks 


p A = E(K S - A)/E(H - K 5 ) 


AJA Ks 


/ 


1.240 


-3.05 1 ± 0.H4 


2.525 ± 0.063 


-3.081 ± 0.011 


2.541 ± 0.006 


H 


1.664 




1.550 




1.550 


K s 


2.164 




1.000 




1.000 


[3.6] 


3.545 


0.649 ± 0.009 


0.643 ± 0.007 


0.558 ± 0.006 


0.69 1 ± 0.004 


[4.5] 


4.442 


0.820 ±0.01 1 


0.549 ± 0.007 


0.754 ± 0.008 


0.582 ± 0.005 


[5.8] 


5.675 


0.985 ± 0.013 


0.458 ± 0.010 


0.866 ± 0.008 


0.520 ± 0.005 


[8.0] 


7.760 


1.023 ± 0.017 


0.437 ±0.011 


0.882 ±0.011 


0.5 12 ± 0.008 
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Fig. 4. fii ow and f3 high for B59 as a function of (H - K s )u mi , for each IRAC band, [3.6], [4.5], [5.8], [8.0] from left to right, showing 
no evidence for a break in the extinction law with column density. The scale in the y-axes is the same for all panels. 




oo- 0.5 



Fig. 5. Pi ow and p high for FeSt as a function of (H - K s )i imit for each IRAC band, [3.6], [4.5], [5.8], [8.0] from left to right, showing 
no evidence for a break in the extinction law with column density. The scale in the y-axes is the same for all panels. 



According to the tests with synthetic data from Paper I, a break 
in the extinction law would show as a significant difference be- 
tween the black (low extinction) and the red (high extinction) 
curves, which is not observed. Based on this analysis we con- 
clude that the departure of the datapoints at high extinction from 
the LINES fit observed in the color-color diagrams is not statis- 
tically significant. We do not find any change of the extinction 
law with column density for either core. 



4. Discussion 

4.1. Comparison with previous work: a universal MIR 
extinction law? 

In the near-infrared the extinction law is well described by a 
power-law (A a °c A~ a ) with indexes tightly ranging from 1.6 
to 1.85. This constancy has led authors to use the near-infrared 
as anchor to normalize the extinction law in the MIR. Although 
early studies reported that the MIR extinction law continued as a 



power-law from the NIR ( |Whittet et al.|1988[|Martin & Whittet 
1 1990 Landini et al. 1984 1, most recent studies find that it departs 
from the power-law around 3 yum and becomes flat (also referred 
to as gray) at least until the silicate feature around 9.7 fim. 

Figure [6] summarizes the extinction laws determined by pre- 
vious works, and shows that most regions have flat extinction 
laws. Each panel shows Aa/Ak s for each of the IRAC bands 
([3.6], [4.5], [5.8], and [8.0]) as a function of the correspond- 
ing reference in the literature. For the same reference, a gray or 
flat extinction law will show approximately the same value in 
all panels, whereas an extrapolation of the NIR power-law will 
show as a decrease from the left- to the right-hand panels. |Lutz| 
[eTaTI ([1996) and |Lutzl ( [1999] Ref. #1) first found a flat MIR 
extinction law toward the Galactic Center using spectroscopy 
of hydrogen recombination lines. Viehmann et al. (2005) later 
measured a consistently flat extinction at 4.66 ^m toward the 
same region, and iNishiyama et al. (2009 Ref. #17) and Fritz 
|et al.| ( |20ir| Ref. #27) independently confirmed a flat MIR ex- 
tinction law toward that line of sight. Based on Spitzer photom- 
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Fig. 6. The extinction law for FeSt (blue filled squares) and B59 (red filled squares), in comparison with the extinction laws from 
the literature, shown for each IRAC band separately for clarity. The scale in the y-axes is the same for all panels. From left to right, 
AxIAk s for the IRAC bands [3.6], [4.5], [5.8], and [8.0] is shown for each reference in the literature. For the works not based on 
IRAC data the plotted values were the closest match in wavelength. The dark circles refer to extinction laws derived toward lines of 
sight containing dense regions, the light circles toward the diffuse medium, the open circles toward low (column) density lines of 
sight in the dire ction of star formi ng regions, and the star symbols toward regions dominated by outflo ws. Also shown are the A^/Ar 
values from the |Draine| ( |2003a|b"] l models for Ry = 3.1 (solid line) and Ry — 5.5 (dotted line), and the Weingartner & Draine] ([2001) 
models for Ry = 5.5, case B (dashed line) for the corresponding wavelengths. All extinction laws are normalized to AhIAk s = 1.55. 

Th e reference s ar e as follows: 

(l)|Lutz]<fT999l l; (2) |Bertoldi etaT] (1999) and |Bertoldi eTaLl|T999) ; (3-5) |Indebetouw et al.|<[2005> for RCW49; I = 284°; I = 42°; (6-10) Flaherty 
etal.| ( |2007) for 2024/2023; 2068/2071; Serpens; Orion; Ophiuchus; (11-12) McClur e]p009|> for A Ks < 1 m ag; A Ks > 1 mag; (13-16) |ZasowskT 
eTar|pn09| > for 10° < |/| < 15°; 30° < \l\ < 40°; 50° < \l\ < 60° ; \l\ > 90°; (17) |Nishiyama et al.|j2009) ; (18-21) plapman et al.lpTO) for 



<A Ks < 0.5; 0.5 < A Ks < 1; 1 < A K < 2; A Ks > 2 mag; ( 22-25) |Chapman & Mundy|(20()9) for < A K < 0.5; 0.5 < A Ks < 1; A Ks > 1 mag: 
regions with outflows; (26)|Gao et"aT]j2009) ; (27) |Fritz etaT]^2011) ; (28-29) |Qlofsson & Olofsson| ( |2011| > for low and high column densities; 
(30-31) |Cambresy et al.| ( |201 1| > for Ay < 15 mag; A v > 20 mag; (32-33) This work for B59 and FeSt. 



etry, Indebetouw et al. ( |2005 1 found a flat MIR extinction law 
toward other lines of sight crossing both dense (Ref. #3) and dif- 
fuse (Ref. #4-5) media. A number of other studies using Spitzer 
followed, namely Flaherty et al.|(|2007| Ref . #6-10), Chapman 
etaLl ( [2059l Ref. #18-21) and |McClure| ( |2009l Ref. #11-12), who 
found gray extinction laws for star forming regions and molecu- 
lar clouds, and Roman-Zunig a et aL](|2007|), [ Chapman & Mundy 
( [2009] Ref. #22-25) and |Qlofsson & Olofsson| ( |2011| Ref. #28), 
who found gray extinction laws for cloud cores. Specifically tar- 
geting low density lines of sight, Zasowski et al. (|2009 Ref. 
#14-16) found flat MIR extinction laws over a very wide por- 
tion of the Galactic midplane. The extinction laws of FeSt (red, 
Ref. #32) and B59 (blue filled circles, Ref. #33) follow the same 
flat trend in the IRAC bands. The column densities probed by 
our data (A v up to 50 mag) are among the highest ever stud- 
ied in this context, and correspond to high local volume densi- 
ties as attested by the presence of high-density tracer molecules 
(Rath borne et al.|2008[|Frau et al.|2010|[20T2) . Still, despite the 
difference in density, the extinction laws of these cores are unre- 
markable with respect to less dense environments. 



So me authors ( |Chapman et al. (2009| R ef. #18 -21), |McClure 
p009l Ref. #11-12) IQlofsson & Olofsson| ( |20ni Ref. #28-29), 
and Cambresy et al. (201 1 Ref. #30-31)) surveyed lines of sight 
of different column densities within single regions and found a 
slight flattening trend toward higher densities, but only the first 
find significantly steep laws in absolute terms. In contrast, other 
studies have looked for but failed to find any significant varia- 
tion of the extinction law with increasing density (Chapman & 
Mundy||2009l Ref #22-24), as do we for B59 and FeSt (Sect. 



3.3| >, even though we probe a large range in extinction. 

There are only three exceptions in the recent literature where 
the observed MIR extinction laws are steep enough to follow the 



gions of Perseus, Ophiuchus and Serpens; Bertoldi et al. (1999) 
and |Rosenthal et al. ( 2000| Ref. #2), toward th e outflow region 
of the Orion OMC- 1 ; and ]Chapman & Mundy| ( [2009l Ref. #25), 
also in the direction of outflows. These differences are interest- 
ing and require further study but represent a small minority of 
the studies in the literature. 

From the data collected so far, there is no clear indication 
that the MIR extinction law varies notably from region to region, 
or that it varies as a function of density, since diffuse and dense 
regions alike show flat extinction laws. The MIR extinction law 
does not seem to depend (uniquely) on density, at least within 
the (column) density regimes probed by current observations. 

4.2. Comparison with models: implications for grain sizes 

As mentioned briefly in the Introduction, the extinction can 
be used to study the grain size distribution of the interstellar 
dust. The extinction law in the optical and ultraviolet has been 
reported to change from region to region, presumably due to 
the size distribution of the grains causing the extinction (e.g., 
Savage & Mathis|[l9"79"} |Whittet|[200? ), with a con-elation be- 
tween larger grains and higher densitieqjinterpreted as evidence 



for grain growth (Card elli et al.|[1989| l. The extinction law at 
these wavelengths is however only sensitive to grain sizes less 
than ~ 0.2 pm (Dr aine|[201 1[ ). To probe the size distribution of 
larger grains, the infrared extinction law must be used instead. 

The observed flat MIR extinction laws are best described by 
the case B, R v of 5.5 model of Weingartner & Draine (2001 



Case B refers to a grain size distribution that includes a signif- 
icant fraction of large (up to 10 pm) grains, whereas their case 



power-law extrapolation from the NIR: Chapman et al. (2009 



Ref. # 18), toward the low column density (Ak, < 0.5 mag) re- 



5 "Larger" grains and "higher"densities refer here to the larger grains 
that can be probed with optical and UV extinction laws (< 0.2 pm) and 
to the largest column densities that can be probed with optical and UV 
observations (A v of a few magnitudes). 



6 



J. Ascenso et al.: The mid-infrared extinction law in the darkest cores of the Pipe Nebula 



A model includes grains with maximum sizes of ~ 1 //m and 
produces a steep MIR extinction law|^] Therefore, flat extinction 
laws are believed to be produced by dust containing significant 
fractions of large grains. Figure [6] shows that the case B, Ry of 
5.5 model of WD01 adequately describes the observed extinc- 
tion laws at 3.6 yum and 4.5 yum but that it does not produce a 
gray enough extinction law to match most observations at 5.8 
pm and 8.0 fim. Since most of these observations are based on 
broad-band data, this may indicate that there are narrow, unre- 
solved emission features contributing to the extinction law at 
these wavelengths. Alternatively, it may mean that the grain size 
distribution corresponding to these models still does not contain 
enough large grains, or large enough grains to explain the obser- 
vations. 

As mentioned before, grains are believed to grow at the high 
densities and low temperatures of shielded molecular clouds, 
where the formation of ice mantles and grain-grain collisions 
are most favored. However, observations show that the extinc- 
tion law is not (uniquely) dependent on density (see Figure [6] 



and Sect. 4.1 



so either other factors determine the shape of the 
MIR extinction law, or large grains can form and exist at rela- 
tively low densities. Even though the case A, Ry = 3.1 model of 
WD01 (later renormalized by |Draine (2003a b)) without large 
grains has been considered to be the appropriate model for the 
diffuse interstellar medium, it is actually hardly ever observed in 
the MIR. 

The extinction laws of B59 and FeSt are also comparable to, 
and actually flatter than the case B, Ry = 5.5 model of WD01 
(Figure|5]l, suggesting that these cores contain a significant frac- 
tion of large grains if the models are taken at face value. Unlike 
Cambresy et al. (2011 1, who find a flattening of the extinction 
law at a fhreahold of Ay = 20 mag, we do not find any such 



change in the extinction laws of B59 or FeSt (Sect. 3.3 i even 



though our data cover a large enough range of extinction and 
our method has been shown to be able to detect it (Paper I), fur- 
ther suggesting that large grains are present already in the low- 
density edges of these cores. However, this is somewhat incon- 
sistent with the results of Aguti et al. (2007 ), who find depletion 
of C 18 toward FeSt at relatively low extinctions {Ay = 12 mag). 
If the accretion of ice mantles by grains as molecular species 
condense onto their surfaces increases their sizes and promotes 
their growth, the depletion of molecules should be regarded as 



indirect evidence for grain growth. For completeness, Frau et al. 
(2012) suggest that there is no depletion in B59, although their 



arguments are solely based on a high abundance of CSQ 

It is clear from this study and from the literature that we 
still do not have a consistent understanding of the extinction law 
and/or of the dust grain properties in the different environments. 

4.3. Comparison between the two cores 

In detail, the MIR extinction laws of FeSt and B59 are not identi- 
cal: the extinction law of FeSt is grayer than that of B59 at the 5- 
to 12-cr level depending on the wavelength. This is true as long 
as our assumption that the NIR extinction law is the same for the 
two lines of sight applies. When converting the reddening laws 
(the color excess ratios f3,\, Sect. 3. 1 1 to extinction laws (A^/A^), 
we assumed A H /A Ks = 1.55 following Indebetouw et al. (2005 ), 



6 Cases A and B are indistinguishable from each other at optical and 
UV wavelengths. 

7 Unpublished radio data of B59 suggest that there may be depletion 
of C 18 at column densities higher than A v = 30 mag (Roman-Zuniga, 
private communication). 



which corresponds to a NIR extinction law being a power law 
with index 1.65 (A^ oc /U L65 ). This value is adopted by many 
studies of the extinction law as universal but in fact different 
values have been measured toward different lines of sight. The 
values in the literature range between 1.85 and 1.65 (see refer- 
ences in the Introduction), and there is no understanding of how 
it changes from region to region. It is therefore possible that B59 
and FeSt have different extinction laws in the NIR. In particular, 
if the index of the NIR extinction law for FeSt is the adopted 
1.65, but that for B59 is instead 1.90, the two extinction laws 
are indistinguishable within the errors. However, as mentioned 



in Sect. 3.2 we measure very similar A j/Ak s values for the two 



cores, supporting the assumption of similar NIR extinction laws 
for the two cores. In the following we take this assumption as 
correct and discuss the implications of the MIR extinction laws 
for the two cores being different, as shown in Figure [3] 

We also assume, for the sake of argument and based on the 
models of WD01, [Drainel ( |2003a| ) and |Draine| ( |2003b| ), that the 
flatter the extinction law, the larger the grains causing the extinc- 
tion. In this context, FeSt should have larger grains (or a larger 
fraction of large grains) than B59, since it has the flattest ex- 
tinction law. The depletion of C ls O at relatively low extinctions 
observed in FeSt in contrast to the lack of depletion in B59 fur- 
ther suports the presence of larger grains in FeSt. 

Why are the grains smaller in B59 than in FeSt? Since the 
two cores belong to the same cloud it is reasonable to assume 
that they have the same chemical composition, so we do not ex- 
pect that different metallicities are the source of the different ex- 
tinction laws. Our data probe approximately the same range of 
column densities for both lines of sight, and both FeSt and B59 
are dense cores, so the column density is likely to represent ap- 
proximately the same physical density. If this is the case, then 
the difference in grain size distribution for the two cores cannot 
be explained by density arguments either. If density is not the 
dominant factor, it may be that the different internal conditions 
of FeSt and B59 are causing the difference in extinction laws. 
In fact, the most obvious difference between the two cores is 
that FeSt is quiescent, whereas B59 contains ~ 20 young stellar 



objec ts ( |Brooke et aT1|2007[ |Forbrich etall2009| |Covey et al 
2010|) and at least two molecular outflows (fOnishi et al.| |1999 ; 
Riaz et al.|2009[ Roman-Zuniga et al. 2010|[Duarte-Cabral et al 



2012). The young stars necessarily increase the temperature of 



the core, at least locally, and their winds and outflows inject me- 
chanical energy and turbulence into the intra-core medium. It 
would therefore not be unexpected that the newly formed stars 
have some bearing in the reprocessing of the dust grains. If star 
formation is instrumental in the destruction of large grains, then 
the extinction law in B59 would be steeper than that of FeSt, as 
is observed. This would be consistent with the other observation 
of steep extinction laws toward other, more exposed outflow re- 
gions (|Bertoldi et al.||1999[ [Rosenthal et al.||2000| |Chapman & 
|Mundy|2009] T~ 



We note that |Alves et al.| ( [2~008| ) find a magnetic field strength 
toward FeSt around 65 /iG, almost 4 times higher than that to- 
ward B59, corresponding to degrees of polarization reaching 
15% in the region of FeSt and around 1% for B59. However, 
the strength of the magnetic field is related to the alignment of 
the grains, and therefore is not very revealing to the understand- 
ing of the grain size distribution. Also, since this study was done 
with R-band observations, it characterizes the magnetic field of 
the lower-density regions (Ay < 5 mag) of the cloud, so the am- 
bient magnetic field instead of that of the cores 
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5. Conclusions 

We applied the unbiased LINES method ( |Ascenso et al.|[20T2] > 
to determine the extinction law from photometric data from 1 .2 
to 8.0 yum toward the two densest cores in the Pipe Nebula, B59 
and FeSt 1-457. 

We find a flat or gray extinction law for both cores. The com- 
parison with models suggests that the dust in these cores contains 
a significant fraction of large grains. Despite the wide range of 
column density probed in this study (up to 50 magnitudes of 
visual extinction) with respect to most previous studies, the ex- 
tinction laws we derive are not remarkably different from others 
toward other lines of sight. In particular a review of the litera- 
ture indicates that the mid-infrared extinction law is gray already 
at low column densities, suggesting that either large grains are 
present already in low column density regions and in the diffuse 
medium, or that the existing dust models need to be revised at 
mid-infrared wavelengths. The implications for the theories of 
grain growth, that generally require high densities for the forma- 
tion of large grains, are not discussed in this paper but should be 
addressed appropriately in the future. 

Although the extinction laws are gray for the two cores, that 
of FeSt is grayer than that of B59. Taken at face value this sug- 
gests that FeSt has larger grains, or a larger fraction of large 
grains, than does B59. We hypothesize that this reflect the inter- 
nal differences in the two cores, in particular with the presence 
of young stars and outflows in B59 and lack thereof in FeSt. We 
propose that the grains have already started to be reprocessed in 
B59 by the increase in temperature and turbulence caused by the 
young stars. 

We find no evidence for a flattening of the extinction law as 
a function of column density, and conclude that the low-density 
(Ay < 10 mag) edges of the cores already contain large grains. 
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